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We report a systematic measurement of the space charge effect observed in the few-ps laser pulse regime 
in laser-based solid-state photoemission spectroscopy experiments. The broadening and the shift of a gold 
Fermi edge as a function of spot size, laser power, and emission angle are characterized for pulse lengths of 6 
ps and 6 eV photon energy. The results are used as a benchmark for an TV-body numerical simulation and 
are compared to different regimes used in photoemission spectroscopy. These results provide an important 
reference for the design of time- and angle-resolved photoemission spectroscopy setups and next-generation 
light sources. 
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I. INTRODUCTION 

Electrons traveling in vacuum to a measurement ap- 
paratus will repel each other because of their mutual 
Coulomb interaction. If a cloud of electrons is sufficiently 
dense, this vacuum space charge effect can alter the time, 
energy, and spatial spread of the cloud, limiting exper- 
imental resolution and resulting in systematic measure- 
ment errors 

The space charge effect has been addressed mostly for 
the design of electron guns for time-resolved electron 
diffraction experiments 2 -^ and has been discussed in 
terms of iV-body numerical simulations, mean-field mod- 
els, and fluid dynamics. Recently the space charge effect 
has also been discussed in angle-resolved photoemission 
spectroscopy (ARPES) experiments 7 ^ 9 , where it has been 
shown to depend primarily on four parameters: photon 
energy, pulse duration, incident beam spot size, and the 
number of electrons photo-emitted per pulse. These pre- 
vious studies explored two different regimes: i) the rel- 
atively high photon energy (> 35 eV) and long pulse 
duration (« 60 ps) 9 regime, and ii) the very low pho- 
ton energy (two-photon photoemission with 3 eV pho- 
tons) and short pulse duration (« 40 fs) regime^ The 
shape of the electron cloud is very different in these two 
regimes. In the first case, the wide distribution of elec- 
tron kinetic energy and long pulse duration produce a 
cloud extended mainly along the direction normal to the 
sample surface, while in the second case the narrow elec- 
tron kinetic energy distribution and short pulse duration 
produce a disk-shaped electron cloud parallel to the sam- 
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pie surface. The former case of a more three-dimensional 
cloud requires numerical simulations while a disk-shaped 
cloud can be characterized analytically using a simple 
power law*£ — 

With the recent advent of laser-based photoemission 
experiments, the need to extend these studies to an in- 
termediate regime becomes imperative. A recent study 
of the space charge effect using laser based ARPES with 
7 eV photon energy and 10 ps long pulse suggests that 
the space charge effect plays a minor role in this regime 
thanks to the mirror charge effect that nearly cancels the 
space charge effect in the 1 to 10 ps pulse duration J£ 

Here we expand these results by reporting a systematic 
study of the space charge effect on polycrystalline gold 
in laser-based photoemission spectroscopy in the inter- 
mediate regime of moderate pulse lengths (« 6 ps) and 
low photon energy (6 eV). We find in the limit of small 
spot size (~ 20 /im diameter) a shift of the Fermi edge 
on the order of 2 meV and a broadening larger than 10 
meV starting already at 100 electrons per pulse. More 
generally, we find that a single power law cannot describe 
the broadening or the Fermi edge shift as a function of 
electrons per pulse and spot size in the various ARPES 
regimes. To explore the non- trivial relation between the 
various experimental parameters, we modeled our data 
with a molecular-dynamics model which employs full TV- 
body numerical simulations of the electron propagation 
from the surface to the detector £ This model is able to re- 
produce quantitatively the space charge effect seen in the 
long and short pulse regimes as well as in the intermedi- 
ate regime explored here. The large range of applicability 
of this model allows one to identify the functional depen- 
dence of the spectral shift and broadening on the number 
of electrons per pulse and spot size. We found that one 
of the key parameters affecting the different power laws 
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^^J_-jU^|J ^ 54.3/n MHz, 150 fs, 840nm 
BBO 2mm 

FIG. 1. (Color online) Laser- ARPES experimental setup. 
The focal lengths are: Ll=10cm, L2=20cm and L3=30cm. 



observed in different ARPES regimes is the maximum 
kinetic energy of the initial photo-electron distribution. 
We believe that this study will be crucial for guiding next 
generation high resolution setups. 



II. EXPERIMENTAL TECHNIQUE 

Experiments were performed on a polycrystalline gold 
sample scratched in air and inserted through a load lock 
into a UHV analysis chamber at low temperature (15-20 
K). Pressure remained below 5 x 10 -11 Torr throughout 
the experiment. In Fig. 1 we show a diagram of the 
experimental setup. The sample geometry is detailed in 
the inset. The spot size was characterized by using a 
wire mesh mounted on the same sample holder as the 
gold sample. 

The laser system is a cavity-dumped, mode-locked 
TiiSapphire oscillator (Coherent Mira), pumped with a 
6 W frequency-doubled NdiYVCU laser. The oscillator 
generates « 150 fs pulses at 840 nm at a repetition rate 
of 54.3/n MHz (n = 1,2,3,...) with pulse energies as 
high as 100 nJ. The photon energy of the laser beam 
was frequency-quadrupled via cascaded, type-I phase- 
matched second harmonic generation in two beta barium 
borate (BBO) crystals of 2 mm and 5 mm thickness. The 
resulting 210 nm (5.9 eV) pulse was estimated to be « 
6 ps long due mainly to group velocity mismatch. The 
UHV chamber window is made of UV- grade fused silica. 
We adjusted the pulse energy by tuning the power ap- 
plied to the cavity dumper Bragg cell. The focusing lens 
(L3 in Fig. 1) has a 300 mm focal length. We varied the 
spot size by defocussing L3. 
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FIG. 2. (Color online) Spot size measurements, (a) Photo- 
electron yield as a function of laser power, (b) Typical x and 
y scan of the integrated photoemission intensity from the W 
mesh. The y scan is vertically shifted 



Assuming a work function of about 4.4 eV, electrons 
at the Fermi energy exit the sample with 1.5 eV kinetic 
energy, corresponding to a velocity of 0.7 /xm/ps. With 
a pulse lasting 6 ps, the initial electron spread normal 
to the surface is thus about 4 /xm. With a spot size 
diameter of about 20-100 /xm, the initial electron cloud is 
a fiat cylinder. This is an intermediate case between the 
fs laser-based data and the synchrotron-based data. 8,9 

The integrated photoemission signal was measured us- 
ing a Phoibos 150 mm hemispherical analyzer from Specs 
with a 0.2 mm entrance slit and 4 eV pass energy, yield- 
ing an energy resolution of 4 meV. The sample was elec- 
trically grounded through a picoammeter, which allowed 
the measurement of the total photo-electron yield. All 
the spectra were averaged over ^18° along the analyzer 
slit except where explicitly stated otherwise. The angle 
between the analyzer and the laser beam is fixed at 55° 
and the light incident on the sample is s-polarized. 

Fig. 2 shows the characterization of the pulse energy 
and spot size. The total photo-electron yield as a func- 
tion of the laser power is shown in Fig. 2(a). The repe- 
tition rate was set to 5.4 MHz. The solid line is a linear 
fit giving a yield of 1 x 10 -4 e~ /photon. For compari- 
son, the typical yield for a photon energy of 35 eV mea- 
sured at BL 10.0.1 at the Advanced Light Source (ALS) 
is about 0.2 e~/photon (1 x 10 12 photons/s with w 300 
nA yield). 9 The difference can be mainly attributed to 
the range of electronic states probed: 6 eV photons only 
excite states in the Au 5d-6s valence band with binding 
energy less than about 1.5 eV, whereas 35 eV photons can 
excite the entire valence band. 9 Additionally, in the lat- 
ter case secondary electrons amplify the photo-electron 
yield noticeably* 1 ^ 

The spot size was measured by scanning a 100x100 
lines/inch tungsten mesh at the sample position while 
measuring the integrated photoemission intensity. The 
mesh wire diameter was 25 jam. The angle (see Fig. 1) 
was set to 55° during the mesh scan so that the x and y 
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FIG. 3. (Color online) Fermi edge of the Au sample as a 
function of the number of photo-electrons per pulse (TV). The 
spectra are normalized and shifted vertically. The solid lines 
are fits. 



scanning directions were perpendicular to the laser beam. 
Fig. 2(b) shows a typical scan along the x and y direc- 
tion obtained from this procedure. We then estimated 
the spot size by deconvolving the peak's full width at 
half maximum (FWHM) using the known mesh wire di- 
ameter. The reported diameter is the square root of the 
x and y FWHM product, corresponding to the diameter 
of a disc of equal area. The ratio between the FWHM 
along x and y is not constant but was smaller than 3. 
We note that with an incident Gaussian beam of 2 mm 
diameter, the diffraction limited spot size with a 300 mm 
focal length lens is about 15 /im. For the FWHM along 
the x direction, we used a corrected wire diameter to ac- 
count for the fact that only the part of the wire facing 
the analyzer could contribute to the measured signal. It 
can be shown that the effective wire diameter in this case 
is 25 x (1 + cos(55°))/2 = 20 /im. The small difference 
in peak separation between the x and y scan is due to 
mesh spacing inhomogeneity. 

Fig. 3 shows typical integrated photoemission spectra 
from the gold sample at various laser powers showing the 
Fermi level broadening and shift. The spot size diameter 
was 18 /im. The edges of these spectra have been fit to 
a function of the form: 

f(E) = a + b- min(£ - E , 0) + c • erfc j 1 '^ ) . 

(1) 

The parameters a, 6, c, Eq, and FWHM are fitting param- 
eters and E is the electron kinetic energy. The first two 
terms account for the background, which includes an off- 



set and a linear contribution for E < Eq. The third term 
is the complementary error function erfc(x), 

2 f x 2 

erfc(x) = 1 - erf (a) = 1 -= / ' dt, (2) 

V 71 " Jo 

which is commonly used to fit Fermi-Dirac distributions 
convolved with experimental resolution. The characteris- 
tic FWHM is defined to be the full-width-half-maximum 
value of the associated normalized Gaussian. 

The top of the panel shows a sharp Fermi edge at low 
power with a FWHM of 7.3 meV, consistent with the ex- 
perimental energy resolution (4 meV from the analyzer 
and 4-5 meV from the laser bandwidth) and the sam- 
ple temperature (about 15 K). The total photo-electron 
yield measurement indicated an average of 1 electrons per 
pulse. As the laser power is increased, a greater number 
of photo-electrons (TV) are emitted per pulse. This results 
in a stronger space charge effect, and the edge becomes 
broader and shifts to higher kinetic energy. 

The shift to higher kinetic energy is due to the fact that 
the electrons at the Fermi energy have the highest ini- 
tial kinetic energy, and therefore the Coulomb repulsion 
of slower photo-electrons can only push them towards 
higher kinetic energy (neglecting any mirror charge ef- 
fects). We note that within the laser power range used 
here, the erfc(x) fit function captures the experimental 
line shape very well. This was not the case at higher 
laser power (not shown). 



III. EXPERIMENTAL RESULTS 

In Fig. 4 we report the shift of the Fermi edge Aef, 
and its broadening T, as a function of number of electrons 
per pulse (TV) for several values of the spot size diame- 
ter (D). Here A^ is defined by the difference between 
the measured edge position (Eq) and the Fermi energy 
(Ep) while T is defined as the edge FWHM deconvolved 
from the sample temperature and experimental resolu- 
tion. The T error bars are mainly estimated from fit 
standard deviation and temperature measurement uncer- 
tainty (estimated to be about 3 K). The Aef error bars 
are the standard deviations plus 0.5 meV. The error on 
TV was estimated to be ±(30 + 10%). 

We observe that T and A#^ increase with increasing 
TV and that, as expected, this effect becomes more severe 
as the spot size decreases. The data from Ref. ^ taken 
in a different regime (hv = 35 eV and pulse length = 60 
ps), corresponding to the largest D value shown in the 
figure (see diamond markers), fall in line with this trend. 
To determine the relation between T, Aef, and TV we 
fit the data with the smallest D value using the same 
power law as the one reported in previous studies^ It is 
clear that neither the linear (see dashed line) nor the \/~N 
dependence (dotted line) can describe effectively the data 
here reported in the ps regime. To model the transition 
from the linear to the \fN regime we fit the data of Fig. 
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FIG. 4. (Color online) Space charge effect as a function of the 
number of electrons per pulse (N) for various spot diameters 
(D). (a) Broadening (T) of the Fermi edge as a function of D 
and N. The dashed line shows a linear fit and the dotted line 
a V^/V fit to the triangular markers, (b) Shift of the Fermi 
edge (Aef) as a function of N for various values of D. The 
dashed lines show a linear fit of the triangular markers. The 
solid lines are fits to a power law. The data set with diamond 
shaped markers was obtained at the ALS with 35eV photon 
energy and was adapted from Ref. 0. The y-axis error bars 
include the errors on the fit and the instrumental resolutions. 
The x-axis error bar is estimated from the sample current 
measurement resolution. The spectra are averaged over ^18° 
along the analyzer slit with — 55° . 



4 with generalized exponents: 

T(N)=a r x A/*, 



(N) 



(3) 



(4) 



In Figs. 5(a) and (b), we report the dependence of the 
fit coefficients a and /3 on D, respectively. Remarkably, 
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FIG. 5. (Color online) Plot of the the fitting coefficients a 
and P as a function of D. The solid lines in panel (a) show 
the average value of f3v and (3 aef- The solid and dashed lines 
in panel (b) are fits to a power law. 



f3 is essentially independent of D for both Y and Aef- 
More specifically, for the broadening we found fir — 0.7, 
which is an intermediate value between the value (1) re- 
ported for synchrotron pulses 9 and the value (0.5) re- 
ported for fs pulses. 8 For the Fermi edge shift we found 
Paef = 1.1. In both cases the prefactor a increases 
with decreasing D. To characterize the space charge ef- 
fect dependence on D more quantitatively, we fit a with 
a function of the form 



a(D) = 



(5) 



where C and A are fitting parameters. In a log- log plot, 
a(D) is a straight line with slope A. By substituting 
a(D) into Eq. ([3]), we find that the space charge effect 
dependence on the number of photo-electrons N and spot 
diameter D is of the from: 



f(N,D)=Cx — . 
More specifically we found for the Fermi edge shift: 
A EF (N,D) = 0.5 x 



D 11 ' 



(6) 



(7) 



5 



and for the Fermi edge broadening: 

N 0.7 

T(N,D) = Ux^, (8) 

where as in all following formulas Aef(N, D) and 
r(7V, D) are in meV and D is in jiva. 
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FIG. 6. (a) r and Aef as a function of the emission angle 
(j>. (Aef) — 77 meV and (r) = 79 meV are estimated from 
the fit of the angle- integrated spectrum, (b) Photoemission 
intensity as a function of (f). The solid lines are fits with the 
function h((j)) = const x cos 2 (0). 

To quantify the angular dependence, the angle-resolved 
results are shown in Fig. 6. The emission angle ((/>) is 
along the same axis as 6 with = 0° corresponding to 
normal emission. The angular acceptance of the analyzer 
was ~ 25° along the entrance slit. Panel (a) shows the 
strong dependence of Aef and T from <p. We show in 
panel (b) the photoemission intensity measured at « 200 
meV below the Fermi energy. As a guide for the eye, we 
fit Aef^)-) ^(0) and the photoemission intensity with 
empirical functions of the form h((j)) = const x cos 2 <p. An 
offset is added in the case of the photoemission intensity. 

The strong correlation between the Fermi edge broad- 
ening and the photoemission intensity suggests that the 
angular dependence of the space charge effect is mainly 
due to the photo-electron emission angular distribution. 
The effect of the angular integration over the analyzer 
acceptance angle is different for Aef and Y. While the 
shift of the angle-integrated spectrum ({Aef) = 77 meV) 
is close to the average value of Aef (0) over 0, the broad- 
ening of the angle-integrated spectrum ((T) = 79 meV) is 
significantly larger than the average value of T ((/>). This 



shows that the broadening due to the angular dependence 
of Aef(^) can be larger than just the kinetic broadening. 
This angular distribution of the space charge effect has a 
direct effect on the data reported here since the spectra 
were integrated over w 18° with 6=55°. 

Since we reported the space charge effects as a function 
of N, measurements at different emission angles will yield 
different values for the coefficient a. This should result in 
a different value of C. The exact contribution of the angle 
integration window and emission angle on the exponents 
f3 and A is more difficult to estimate. 



IV. SIMULATION RESULTS 

The numerical simulations are based on a representa- 
tive electron energy distribution with a sharp Fermi edge 
at 1.56 eV. The pulse duration was assumed to be rect- 
angular with a duration of 6 ps. The simulations were 
repeated at least 10 times for each parameters setting. 
Aef and T are extracted with the same edge fit func- 
tion with the exception of the parameters a and b being 
set to 0. The electron emission was assumed to follow a 
cos 2 ((f)) distribution following the experimental results of 
Fig. 6 and only electrons with an emission angle between 
46° and 64° are integrated in order to mimic the experi- 
mental condition. The number of electrons per pulse TV 
varied between 200 and 4000. Spot diameters between 
18 jam and 359 jam were used. 

Fig. 7 shows the simulated broadening T and shift Aef 
of the Fermi edge as a function of the number of electrons 
per pulse N for different spot diameters D. Fig. 8 shows 
T and Aef as a function of the spot diameter D for 
different numbers of electrons per pulse TV. We fit the 
results for T and Aef with the same power law as in the 
experimental case and show the comparative results in 
Table 1 . The exponent values are the average value of all 
the fitting results (not shown). 





Simulation 


Experiment 


Aef(N, D) = 


TV - 79 

3-lX £>o.85 


0.5x DlA 


T(N, D) = 


^0.83 
O.OX 2)1-02 


N 0.7 

13x jj 12 



TABLE I. Summary of the measured and simulated exponent 
values. 

The broadening and the shift of the Fermi edge show 
a good overall agreement between simulation and exper- 
iment. Small systematic deviations are observed in the 
small effect regime. In particular for N < 600 and/or 
D > 200 /xm. 
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FIG. 7. (Color online) Simulation of the space charge effect, 
(a) r and (b) Aef as a function of N for various D value. The 
error bars result from multiple simulation runs and represent 
the standard deviations. The solid lines are fits to a power 
law. 
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FIG. 8. (Color online) Simulation of the space charge effect, 
(a) r and (b) Aef as a function of D for various N values. 
The error bars result from multiple simulation runs and rep- 
resent the standard deviations. The solid lines are fits to a 
power law. 



V. DISCUSSION 

In the case of the disk-shaped photo-electron cloud pro- 
duced by a fs pulse, the broadening due to the space 
charge effect was proposed to be of the form 8 




On the other hand, for the 60 ps synchrotron pulse, 
the broadening follows a linear relation with N. The 
T oc TV 7 (iV 8 for the simulation) dependence we re- 
port for this intermediate regime of laser-based ARPES 
with 6 ps pulses is between these two values. However, 
these three regimes cannot be solely compared on a pulse 
duration basis. The fs pulse analysis is based on a sur- 
face state peak (86 meV FWHM) measured at the band 
minimum in the Brillouin zone center. In this case, the 



majority of the photo-electrons have a higher kinetic en- 
ergy. This will affect the space charge effect in two ways: 

i) the peak position will be shifted to lower kinetic energy; 

ii) the momentum resolution will decrease resulting in a 
broader spectrum shifted towards higher kinetic energy. 
While the total shift in energy ends up being negligible, 
the contribution of the momentum broadening on T is 
not£ This makes a direct comparison difficult but fur- 
ther comparative analysis can be done with our model. 

According to the simulation results, the nonlinear ity 
T oc 7V^ r with /?r ^ 1 seems to depend largely on the 
initial photo-electron distribution. More specifically, the 
simulations show that the exponent f3r decreases as the 
maximum electron kinetic energy decreases. This expo- 
nent is around 1 for high kinetic energies, decreases to 
Py = 0.7 — 0.8 in the intermediate regime presented here 
with a maximum electron kinetic energy of about 1.5 eV, 
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and decreases down to f3r = 0.5 for the fs pulse exper- 
iment where the maximum electron kinetic energy was 
around 1.1 eV*£ We found also that the simulation shows 
no indications that the pulse durations affect this behav- 
ior of the exponent fir- 

Similarly, we find the broadening dependence on spot 
size to be proportional to D~ 12 for the simula- 

tion), compared to D~ 0,5 the fs pulse data, though more 
experimental data are required here to validate the trend. 

It is interesting to note also that the ratio of the two 
coefficients /3r and Ar is not 1 as for the fs pulse case, but 
closer to 0.6-0.8. Since the error bars on A are relatively 
large, fits with X = f3 are shown with a dashed line in 
Fig. 5 (b). From this fit it is clear that additional data 
points are needed to conclude on the correct value of 
the ratio between A and f3. Nevertheless, these results 
set already a limit in the range in which the analytical 
model proposed in Ref. 8 is valid. 

VI. CONCLUSIONS 

In conclusion, we report the measurement of the space 
charge effect in the few-ps laser pulse regime. We ob- 
serve a broadening of the spectral feature exceeding 10 
meV and an energy shift exceeding 2 meV with pulses of 
as few as a 100 photo-electrons for small spot sizes (« 
20 /im diameter). For larger spot sizes (~ 165 /im diam- 
eter) similar effects are observed only as the number of 
photo-electrons per pulse reaches 600 to 800. We found 
that neither the linear model that fits synchrotron data 
nor an y/~N model that fits fs pulses describe our data ef- 
fectively, and we introduced a simple phenomenological 
model with generalized exponents that characterizes the 
broadening and energy shift of the spectra as a function 
of spot size and number of photo-electrons over the en- 



tire range probed. The results were used as a benchmark 
for an TV-body numerical simulation with a good over- 
all agreement. In turn the simulation allowed us to ex- 
plain the difference observed in different ARPES regimes 
in terms of initial electron kinetic energy. These results 
represent a reference for the design and interpretation 
of photoemission experiments with next-generation light 
sources, such as free-electron lasers and table- top laser 
sources. 
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